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Abstract 

Inhomogeneous recombination can give rise to perturbations in the electron number density 
which can be a factor of five larger than the perturbations in baryon density. We do a thorough 
analysis of the second order anisotropies generated in the cosmic microwave background (CMB) 
due to perturbations in the electron number density. We show that solving the second order Boltz- 
mann equation for photons is equivalent to solving the first + second order Boltzmann equations 
and then taking the second order part of the solution. We find the approximate solution to the 
photon Boltzmann hierarchy in t modes and show that the contributions from inhomogeneous 
recombination to the second order monopole, dipole and quadrupole are numerically small. We 
also point out that perturbing the electron number density in the first order tight coupling and 
damping solutions for the monopole, dipole and quadrupole is not equivalent to solving the second 
order Boltzmann equations for inhomogeneous recombination. Finally we confirm our result in a 
previous paper that inhomogeneous recombination gives rise to a local type non-Gaussianity pa- 
rameter /tvl ~ — !• The signal to noise for the detection of the temperature bispectrum generated 
by inhomogeneous recombination is ~ 1 for an ideal full sky experiment measuring modes up to 

^■max — 2500. 
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I. INTRODUCTION 



The process of recombination depends on the energy density of photons and baryons 
as well as the number density of electrons. Perturbations in energy and number density 
of photons, baryons and electrons therefore makes recombination a function of position. 
The resulting perturbations in the electron number density, 5 e , give rise to second order 
perturbations in the photons through Compton scattering. The perturbations in the electron 
number density were first calculated by Novosyadlyj jlj], who found that 5 e ~ 5 x <5 b on large 
scales, where 5b is the perturbation in the baryon density. Recently Senatore et al. [2] did 
a more rigorous analysis, including perturbations in the escape probability of Lya photons, 
and found a similar result. 

The factor of five enhancement of the electron number perturbation suggests the possi- 
bility of observable non-Gaussianity even if the initial conditions are completely Gaussian. 
Assessing whether these effects are observable by Planck 3| is therefore important, espe- 
cially since Planck aims to probe the non-Gaussianities in the initial conditions. There have 



been many studies of different second order effects [4|-2lj]. In our previous paper 22] (here- 
after KW09) we calculated the bispectrum arising due to inhomogeneous recombination and 
found that it gives rise to a local type non-Gaussianity with the non linear (NL) parameter 
\Inl\ ~ 1- However we ignored the second order photon monopole and quadrupole and 
electron velocity in the second order Boltzmann equation. In this paper we justify ignoring 
these terms. We also examine two different methods of arriving at the second order solutions 
to the photon Boltzmann equation. The first method is to solve the first and second order 
Boltzmann equations together and take the second order part of the resulting solution as the 
solution to the second order Boltzmann equation. The second method is to solve the sec- 
ond order Boltzmann equation separately. In KW09 we solved the second order Boltzmann 
equation separately and found that the first order photon monopole does not contribute to 
the second order anisotropy while the first order photon dipole is partially cancelled by the 
first order electron velocity. We prove that the two methods are equivalent. This is also 
important for the self-consistency of the perturbation theory. The important fact that the 
first order source terms are suppressed is somewhat obscured in the expression resulting from 
solving the first and second order equations together. We also explain in the conclusions 
section that perturbing the number density of electrons in the first order tight coupling and 



damping solutions for the monopole, dipole and quadrupole is not equivalent to solving the 
second order Boltzmann equation for inhomogeneous recombination. The method of per- 
turbing the first order solutions was followed in 23J whereas what we want is the solution to 
the second order Boltzmann equation which we find in this paper. Following cosmological 
parameters are used for numerical calculations: baryon density flf, = 0.0418, cold dark mat- 
ter density Q c = 0.1965, cosmological constant Q\ = 0.7617, number of massless neutrinos 
N v = 3.04, Hubble constant H = 73, CMB temperature Tqmb = 2.725, primordial Helium 
fraction t/H e = 0.24, spectral index of the primordial power spectrum n s = 1.0, and ag = 0.8. 
All first order quantities are in conformal Newtonian gauge and calculated using CMBFAST 
24J]. Electron number density perturbation is calculated using DRECFAST [l]. 



II. LINE OF SIGHT INTEGRATION AT SECOND ORDER: METHOD 1 



We begin with the first + second order equations as given in, for example, 25( Equations 
6.6 and 6.11. We drop the second order metric perturbations and products of first order 
terms which do not contain 5 e , the electron number density perturbation. However we 
retain the full first order equation since it gives rise to second order terms, as we will later 
see. We drop the usual factors of 1/2 multiplied with the second order variables, and 
use 0w = Aw/4 as our perturbation variable for convenience. O = 5T/T is the photon 
temperature perturbation while A is the perturbation in the photon distribution function 
integrated over momentum and normalized appropriately 25j. Superscripts (i) denote the 
order of perturbation. In what follows all perturbation variables are functions of coordinates 
on spatial hypersurface x, line of sight angle n and conformal time rj in real space and 
functions of Fourier mode k, n and rj in Fourier space unless specified otherwise. We will 
use same symbols for real space and Fourier space quantities but that should not cause any 
confusion as only one quantity is needed at a time. Boldface quantities are 3-vectors while 
"indicates a unit 3-vector. We use the following metric signature with (f> = </>W + cj)( 2 > + ... 
etc. and ignoring vector and tensor modes 

ds 2 = a 2 (r]) [-e^dr] 2 + e^dx 2 ] (1) 

Also we decompose the first order temperature perturbation in Fourier space into I modes 
as 0^(77, k, n) = J2e(~ i e )(2£ + l)^(fi-k)©^ 1 (rj, k), where P^(n.k) are the Legendre poly- 
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nomials. For the second order temperature perturbation we use the spherical harmonic de- 
composition defined by, 0^(77, x) = / dnQ^jri, x, n)>X,(n) and similarly in Fourier space. 



Also the 



Lit 



Note that this differs from the convention used in [25( by a factor of 
electron velocity, v e ^, is equal to the baryon velocity to a high precision and we will drop 
the subscript on v e in the rest of the paper. 

We start with the first + second order Boltzmann equation for photons in real space, 
ignoring second order metric perturbations and second order terms which are products of 
first order terms but do not contain S e = (n e — n e )/n e , where n e (f], x) is the electron number 
density and n e (rj) is the mean electron number density. 

± [0« (77, x, n) + ^ (r/, x) + 0< 2 > (77, x, n)] - J- (0« (77, x) + ^ (77, x)) = 
n^ara^) \(l + ^(77, x)) (C (1) (r/, x, n) - 0«(r7, x, n) - 0^(r/,x,n 



+ 



' eg^x) + ^£ 0^(77, x)Y 2m (n) + v( 2 >fo,x).n 



(2) 



where we have defined which is given in Fourier space by 

C«(77,k,n) = ef } (»7,k) - ^0f ) (7 7 ,k)P 2 (k.n) +v( 1 )(7 7 ,k).n, 



(3) 



denotes the total derivative which is equal to ^- + n l -£j along the line of sight to zeroth 
order, n denotes the line of sight direction, <jt is the Thomson scattering cross section. We 



now 


add 


" d 

dr] 


— T 



#07,x,n) = 



drj ^ 



0W + + (2) 1 = P(?7,x,n 



(1 + #)) (c?w + + ^=eg> + 1 £ eSr 2m (n) 



where we have defined f(?7) = — n e 0Ya, with r(r/) = — tdrj. tjq is the conformal time at 
a — 1. Now we use the fact that along the photon geodesic x is a function of rj to write 
Equation @] as 



J? *_ [(©(D + ^(D + 0(59) gC^+^K' 



i?(77,x,n) (5) 



Note that the above equation can only be written if the integrals appearing are evaluated 
along the line of sight and so x ceases to be an independent variable outside the integrals. 



Jx(r,) 



Integrating Equation [5] formally along the line of sight results in 

(@(d + ^(D + @( 2) ) UM = r dv j:° «H 1+Sil) ) ^ [R{v , X) a)], 

•/ o 

= fj dqe* (l + jf* drffS^ | x(r)0 ) [R( v , x, fi)] xW (6) 

In the last line we have assumed that 0° drj'fd \ x (rf) is small compared to unity and ap- 
proximately of same order as which is a good enough assumption once recombination 
starts. 

Taking the second order part of the above equation we get 

e (2) u o) (r/o) = 



+ 



(-r) (C7« + + ^| + 1 E eSF 2m (n) + v^.nj 



rvo 




/ dr/e r 


(-r)j 


/o 





(7) 



If we consider a single observer then we don't have an independent three dimensional space 
variable with respect to which we can Fourier transform this equation. If we consider all 
possible observers then y = x(r/o) spans all space at time r] and we can write x(r/) = 
xo + nrj = y + n(?7 — 770) along the line of sight. Now all quantities in Equation [7] are 
functions of the same variable y and we can take Fourier transform with respect to it. 
The result is (Note that all perturbation variables are Fourier transforms of the respective 
quantities in the rest of this section, we omit the arguments (k) where there is no confusion.) 



e< 2 >(7fc,k,n 
no 



:-r(v)) 



( / |^ (k '' V) (k " k '' V) + ^ (k " k '' V) 



u + 15 E Q ?1(V, k)F 2m (n) + v< a >fa, k).n 



/ ~~~~ 7 i r^v k '^'-")f(? 7 , )^ i) (k',7 7 ') 



(2vr) d Jr, 



— I 

dr] 



« (k - k, 77) + v« (k - k, 77)) - f (77) (c (1) (k - k, 77) + v« (k - k, 77) 



(8) 



where we have used the properties of Fourier transform when the variable getting trans- 
formed is shifted and which gives the phase factors on the right hand side. We could also 



have chosen initial point xo = y' or x(^i) = yi as our integration variable for any fixed rji 
and got the same result. 



III. LINE OF SIGHT INTEGRATION AT SECOND ORDER: METHOD 2 

Another way to do the formal integration of the Boltzmann equation is to move all terms 
containing 8)P and primordial potentials to the right hand side in Equation [2j take Fourier 
transform of the resulting equation and then integrate along the line of sight. This is in fact 



what is done in 



25] and KW09. In that case the solution for is 



Jo 



ik.(x(r?)-x(»?o)). 



-T 



( Z |^^ )(k/) " " @(1)(k " k0 )) 







(2) 
,„, 



V47T 1U m 



(2) 



(9) 



We now integrate by parts in variable r] the term involving 0^. The boundary terms vanish, 
resulting in 



drje 
d 3 k 



ik.(x(77)-x(?7o)) e -T^. 



7|^V)e^-k')) 



J H_ e ~i^{vo) £ ^ e *'-*(^«(k') (e- T e^^Q^\k - k!) 



(2tt) 

d 3 k' 
(2tt) j 



-ik.x(?7i 



(10) 



We now use the first order equation for 0^ to obtain 
d 3 k' 



(2vr 



■_ e -i^( V0 ) j* dr] | j* dr/^fn+fflsW (k', rf) } e^ e ^ k - k ')- x W 



-fC (1) (k - k') - i(k - k').m/> (1) (k - k') + 
d 3 k' 



(27T) 



3 ik.(x(r;)-x(77o)) 



J drj^j dVe ik '- (x(?? ' )_x( ' ?)) T(V)^ 1) (k / ,V)}e~ T 



x 



-f C (1) (k - k') - z(k - k').n^ {1) (k - k 



^ , <90(k - k') 
drj 



(11) 



By doing integration by parts once again on terms containing ip in Equation [TTJ similar 



to what is done in solving the first order Boltzmann equation [26] , and then using the result 
in Equation [HI we obtain Equation [SJ This shows the simple connection between the two 
approaches. 

In KW09 we worked with Equation O In Equation it is readily apparent that there is 
cancellation between the collision term C^ 1 ' and This point is somewhat obscured in 

Equation [8] since the cancellation is now happening between 5 e terms. Nevertheless we have 
shown the exact equivalence of the two approaches and that there is cancellation of first 
order terms which leads to a small value of /jvl even though the electron number density is 
enhanced by a factor of ~ 5. It is also clear from Equation that the term which causes the 
cancellation, #WgC0, has no direct counterpart among the source terms in the first order 
Boltzmann equation. Thus we have to be careful while using analogies with the first order 
Boltzmann equation to estimate the second order solutions. We will return to this point in 
the conclusions section. 



IV. BOLTZMANN HIERARCHY AT SECOND ORDER 



The Boltzmann equation for photons in Fourier space, ignoring all the first order terms 
that do not involve the electron number density perturbation is [25[ 

6 {2) (k, n, rj) + m.k9 (2) (k, n, rj) - f0 (2) (k, n, rj) = S (2) (k, n, rj) 



5( 2 )(k,n^)^-r/^«(k-k'^) 



t" 



+n.kV 1 )(k',77) - -P 2 (k / .fi)n( 1 )(k / ,77) 

B (2) 
u 00 



, T X V) + jo E © w (k, v)Y2m> (n) + v( 2 ) (k, V ).h 
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•£ (-,)'" (2f + l)P < »(n.k')e^(k',,) 



i">2 
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e (2) 



00 



k, V) + i5 E w (k, r,)Y 2m , (n) + (k, V ).h 

m' 



(12) 



where V^ 1 ) is the first order photon velocity. V"W and V^ 2 ), the second order photon 



velocity are defined as follows [251 ]: 



(27T) 3 ' 

V 7 ( 1 )(k') = ^/rfn6( 1 )(k',r ? ,n)n, 
V( 2 )(k,r/) = — / drie (2) (k,77,n)n-4 
dn@^ 2 ^(k, 77, n)n 



d 3 k' 
(2^ 



eJ 1) (k-k',r / )V( 1 )(k',r / ) 



4tt 



(13) 



In the last line we have ignored the second term since it does not contain 5^. We remark 
that this extra term in the above equation partially cancels a term of the form 0q x v in the 
full second order equation. The dot product of photon velocities with line of sight direction 
which appears in the Boltzmann equation is given by 

V 7 ( 1 )(k').n = -ie?\k', V )4nJ2Y* ml (k')Y lm ,(h) 

m' 

Vf (k, rj) A = £ (k, rj)Y lm , (n) . ( 14) 



We choose z axis along k and take the spherical harmonic transform of Equation 
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In above we have defined 

d 3 k' 



eg 5 m 5 m0 + ±e&6 a + VW6 a + stm 



Sv 



y^(k') (^ 1) (k / ,77) + 3zeW(k',r 7 ) 



(15) 



(16) 



and Vj^'Sgi is the spherical harmonic transform of v^.n. All second order quantities are 
functions of (k, 77). Note that different m modes are independent of each other. Now we can 
write down the Boltzmann hierarchy explicitly. 
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(17) 
(18) 
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For £ > 3, 
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lm 
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e (2) + 
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\ (2£ + 3)(2£ + l) 



f+lm 



•J, 



(2tt) 3 
d 3 k' 
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(19) 



We note that the first order monopole does not appear in the above equations. Also the first 
order photon dipole is partially cancelled by the first order electron dipole. Thus only the 
first order quadrupole and higher multipoles contribute to the hierarchy. These first order 
terms are small during recombination and thus we should expect the second order terms due 
to inhomogeneous recombination to be small. This cancellation counteracts the production 
of non-Gaussianity due to enhancement in 5^}\ 

V. APPROXIMATE SOLUTION OF BOLTZMANN HIERARCHY 



To find the approximate solutions we can use the fact that during recombination f >> 
1/rj. Then, as in the case of the first order Boltzmann equation, we can attempt to find 
an approximate solution at different orders in 1/f. In the limit of f >> l/rj, which is true 
during the entire recombination period except at the very end when the visibility also drops 
sharply, we can ignore the I > 3 modes. Also in Equation [18] we can ignore terms with I > 2 
which do not involve f. Equation [18] with these approximations is 



B (2) 
^2m 



lOik jA-w? (2) 10 



9t 



15 



(20) 



Using this in Equation [T7] 
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(21) 



To proceed further we need the momentum equation for baryons [271 ] . Note that we 
ignore the second order metric perturbations and the terms arising from the first order 
oerturbations that do not contain 8^ as we did with the Boltzmann equation for photons 



25] 



drj 



+ - 
R 



(v«(k' )?7 )-V«(k', 



+ (v( 2 )(k,r / )-V( 2 )(k,r ? ) 



r 



+ (v( 2 )(k,r ? )-Vf(k,r / ) 



(22) 



We have defined ratio of mean baryon to mean photon density R = 3p&/4p 7 . Ignoring the 
expansion term above introduces only a small error on small scales (factors of (1 + i?) 1//4 ) 
which is not important here (for example see Chap 8, Exercise 5 in 26[, also 28(). We 
take the dot product of above equation with line of sight direction n and take the spherical 
harmonic transform of the resulting equation. The result is 

dV® f 

_ r m J_ 

dr] R 



(2) 
m 



(23) 



We can expand Equation [231 perturbatively in R/f as in the first order case 2y, |28|. At 
zeroth order in R/f all the source terms (terms which are products of the first order terms) 
vanish. This causes all the intrinsic second order terms to also vanish if we impose Gaussian 
initial conditions. Thus all terms in the hierarchy are of first order or higher in R/f. At 
first order in R/f we have 

Vg* = effi - Sf v (24) 
Using this in Equation [23] we get up to second order in ■? 



T/(2) _ R (2) _ qm , ( G (2) _ , m \ 



f drj 



(25) 



Continuing like this we can obtain the terms at higher orders in ^. Note that in first order 
perturbation theory we need to go to second order in factors of j to get the damping solution. 
However here we are interested in the contribution of 5 e to the second order anisotropies 
which are intrinsically of first order in -? and it suffices to work at first order in visible factors 
of This gives us the leading term in the solution of the second order Boltzmann equation. 
We comment on the solution beyond this approximation in Appendix [Bl At leading order in 
-? the equations simplify a lot and the solution is similar to that of the first order Boltzmann 



equation 



Using Equation [25] in Equation [21] we get (dropping a higher order term from 
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Equation [21 
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The solution to this equation in the limit that the sound speed c s = J 1/3(1 + R) is slowly 
varying is given by 



©00 = Ci sin [kr s (r])} + C 2 cos [kr s (rj)] 



drj' 



(28) 



where we have defined the sound horizon r s (jj) = f$ drj'c s (jj'). With the Gaussian initial 
conditions, the second order part of temperature anisotropy and its derivative are initially 
zero. Thus G\ = G% = 0. Integrating by parts the Ss v term we get, assuming slowly varying 
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Taking derivative with respect to 77 of above equation we get 
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For m = ±1 modes we can directly integrate Equation 



Jo dr? 9(l + J R(r / '))Vl5 52 ^ )+ l + i? ^ 



(31) 
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We can combine Equations [30] and [31] to get 
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The quadrupole is given by ignoring the 1 / (r) term in Equation [20] (at the level of approxi- 
mation we are working). 
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Finally the second order baryon velocity is given by (Equation 
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An important point to note here is that the photon and baryon velocities are not equal. 
In particular the sign of the last term above is different (in addition to a factor of R) . These 



were assumed to be equal in 



23]. 



VI. NUMERICAL RESULTS 

We want to calculate the angular averaged bispectrum due to 0qq , Vff and 0|>m- The 
contribution from 0qo as well as the S$2 terms in Vff to the angular averaged bispectrum 
is exactly zero. This is shown in Appendix |A] The reason that the contribution from 0qo 
vanishes is the absence of first order monopole from the second order Boltzmann equations. 

(2) 

The contribution to 0q O from the first order dipole and quadrupole averages to zero. Same 
is true for the contribution from first order quadrupole terms in Vj£\ 

Thus the only terms which will give non-zero contribution to the angular averaged bis- 
pectrum are 0^ an d Sg v terms in 0^ an d the last term in Equation [3H are same 
as the terms already calculated in KW09 with additional multiplying factors. The integral 
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FIG. 1. sin [k(r s (ji) — r s (rf))] for k = 0.25 as a function of if for different values of 77. All curves 
end at 7/ = 77 
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FIG. 2. 3@f> —iv^ as a function of conformal time 77 for wavenumber k = 0.2, 0.3, 0.4 Mpc 1 . Note 
that it becomes almost monotonically increasing at large 77 when photon free streaming becomes 
important. 



term involving S$ v in Equation [M] can be calculated exactly following the calculation in 
Appendix [A] However there is an easier way to estimate the magnitude of this term. Figure 
[1] shows the function sin \k{r s (rj) — r s (r]'))] at k — 0.25 for different values of 77 as a function 
of vf . In general there will be cancellation due to oscillations in the sin [k(r s {rf) — r s (rj'))] as 
well as Ssv (Figure [2] and |2J). We can get an upper bound for the region after the peak 
of the visibility function when the magnitude of 3B^ — iv^ is monotonically increasing 
by assuming that the last half cycle of the sine contributes without any cancellation and 
Ssviv') ~ Ss v (r]). Thus we arrive at the following approximation (with slowly varying sound 
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FIG. 3. Gg , — iv^\ and | as a function of t] for wavenumber k = 0.001 Mpc 1 . Also 

shown is the visibility function g(rj) = —fe~ T . 



.0000 



0.1 000 



0.01 00 



0.001 



0.0001 



/. , , . . . : 

260 280 300 320 340 360 380 400 
7] (Mpc) 



FIG. 4. Gg 1 ^, |30i — iv^\ and |Gg | as a function of rj for wavenumber k = 0.01 Mpc 1 . The key 
is the same as in Figure [3l 

speed assumption) 
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l + Rfry™'"' 1 ""' (35) 
where the < sign is understood to be with respect to the magnitude of the terms. For 
most values of r\ and k, where we don't have a monotonic 30^ —iv^\ there will be additional 
cancellations due to the oscillations in 30^ — iv^\ Thus the above term will be smaller 
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FIG. 5. @ { \ 3@\ > - iv^ and G^ as a function of r\ for wavenumber k = 0.1 Mpc 1 . Note that 
at small scales 30^ — iv^ becomes comparable to Gq 1 ^, but its contribution to the bispectrum is 
suppressed because it is weighted by the derivative of spherical Bessel function. See also Equation 
10 and Figure 3 in KW09. 
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FIG. 6. \ 39^ - i«W and G^ function of r\ for wavenumber k = 0.2 Mpc 1 . Note that 
at small scales 30^ — iv^ becomes comparable to Gq 1 ^, but its contribution to the bispectrum is 
suppressed because it is weighted by the derivative of spherical Bessel function. See also Equation 
10 and Figure 3 in KW09. 

than or at most of similar magnitude as the last term in Equation [3D As we will see later 
the last term in Equation [34] gives only ~ 5% contribution to signal to noise and is thus not 
important. 

Before presenting the numerical results we note that S$ v remains small until the very 
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end of recombination. By the time Sg v finally becomes somewhat larger the visibility 
function becomes small suppressing the contribution to the CMB anisotropics. Figures 
and show comparison between Gg 1 , 3G^ — iv^ and G^ for wavenumbers k = 
0.001 Mpc -1 , 0.01 Mpc -1 , 0.1 Mpc -1 and 0.2 Mpc" 1 . In interpreting these figures it should 
be kept in mind that 3G] 1 — iv^ is weighted by the derivative of the spherical Bessel func- 
tion (Equation 10 in KW09) in the expression for bispectrum which is smaller than the 
spherical Bessel function by about an order of magnitude near the peak. Thus even though 
in Figure El and El 36 S 1} — iv^ seems comparable in magnitude to its contribution to 
the bispectrum is much smaller. 

We will collectively refer to the source terms calculated in KW09 as S KWog , that is all the 
terms on the right hand side of Equation El except Q^J and Gj>^. Figure [7] shows the 
confusion with primordial bispectrum of local type as parameterized by /jvl defined in KW09 

(2) 

as a function of maximum £ mode measured by an ideal experiment due to Q^m — 10/95^ 

and = -1/(1 + R)SR. For i max = 2500 we get f NL 0.02, a few percent of the value 

found in KW09 for S KWog . An important point to note is that the sign of the bispectrum 
at small scales is same as the net contribution from S KWog . Thus the new terms calculated 
here will add to the bispectrum from S KW09 and should increase S/N by a small amount. 

In Figure [8] we show the signal to noise ratio for the detection of the bispectrum generated 
by inhomogeneous recombination for a cosmic variance limited experiment as a function of 



the maximum multipole moment observed £ m ax 29] 

S 1 



- ■ 



rp _ \ " K^rec ) 



Afo,^ = 1 + S^iz + 6t 2 i 3 + + 28i 1 i 2 8e 2 £ 3 , (36) 

where B^ 1 ^ 3 is the angular averaged bispectrum generated by inhomogeneous recombination 
, Ct is the CMB angular power spectrum and 5e 1 e 2 is the Kronecker delta function. We get 
S/N ~ 1 at £ max = 2500. Contributions from S KWog and an d Vff calculated in this 
paper are shown separately. S KWog give S/N ~ 1 compared with S/N ~ 0.05 contributed 
by the second order baryon velocity and second order quadrupole. A future high resolution 
cosmic variance limited experiment may thus see a hint of inhomogeneous recombination in 
the bispectrum. 
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FIG. 7. Confusion with primordial non-Gaussianity parameterized by /jvx- Contribution of Q^m = 
10/95^ and vJa ] = -l/(l+R)Sf v is only a few per cent of the contribution from S KW09 , the source 
terms calculated in KW09. S KWm giv6S ci cumulative contribution of Jml ^ — 1 &t ^rnax — 

2500. 

The calculations were done including Fourier modes up to k = 0.5 Mpc -1 . Contributions from 
k ^ 0.4 Mpc -1 are negligible. 
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FIG. 8. Signal to noise ratio for the bispectrum generated by inhomogeneous recombination for a 
cosmic variance limited experiment as a function of the in ax i nm m multipolo moment ^rnax- 

S/N due 

to S KWm is ~ 1 for l max = 2500. Contribution due to 8^ = 10/95^ and vJh } = -l/{l+R)Sf v is 
only a few percent of the contributions S KW09 . Also shown for comparison is S/N from primordial 
non-Gaussianity with /jvx = 1- The calculations were done including Fourier modes up to k = 
0.5 Mpc -1 . Contributions from k > 0.4 Mpc -1 are negligible. 
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VII. CONCLUSIONS 



We have analyzed two different ways of integrating the second order photon Boltzmann 
equations. It is necessary for the consistency of perturbation theory that it should not 
matter if you solve different perturbation orders together or separately and we find that it 
is so in this case. We can define a typical second order term to be of the form x 0g 
with a prefactor of order unity and which can be expected to give rise to a local type non- 
Gaussianity parameter |/jvl| ~ 1- Then we have shown that the second order monopole, 
dipole and quadrupole are smaller than typical second order terms. Although we have 
derived this result in the tight coupling limit to second order in R/f, the fact that these 
terms are small is valid in general. This is because the cancellation that causes these terms 
to be small occurs in the original Boltzmann equations. 

It can be seen that perturbing the electron number density in the first order monopole, 
dipole and quadrupole solutions does not work as follows. The full first order solution can be 
approximately written as a product of an oscillating part and a damping part. Senatore et 



al. 23] perturb just the damping part to estimate the second order solution. The oscillating 
part of the solution does not contain explicit dependence on the electron number density but 
;he equations used in arriving at that solution do depend on the electron number density 



261 ] . To get the oscillating part we have to expand the baryon momentum equation to first 
order in R/f. The factor of f however cancels when the baryon momentum equation is 
substituted into the photon Boltzmann equation and does not explicitly show up in the 
resulting oscillating solution. Similar cancellation happens for the damping solution as well. 
When the electron number density is perturbed in the original equations these additional 
factors of f lead to additional terms in the second order equation that depend on electron 
number density perturbation. Thus there is no way to perturb the electron number density 
in the first order oscillating and damping solutions to take into account these extra second 
order terms and the only way to get the correct second order solution is to solve the second 
order Boltzmann hierarchy explicitly as we have done. In particular the terms missed come 
from the 5 e Q^ term in the second order Boltzmann equation which also results in the 
cancellation of the first order monopole in the second order Boltzmann hierarchy and gives 
the second 5 e term in Equation 

In addition the correct solution should satisfy the relation between the second order 
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monopole and dipole, Equation [T71 (first equation in the Boltzmann hierarchy). The solutions 
given in Senatore et al. [23] clearly fail to satisfy this relation. In particular this relation says 
that the second order monopole and dipole should have the same dependence on angular 
wavenumbers, the factors of Y^ m (k). The first order solutions are the solutions for the 
transfer functions and depend on only the wavenumber magnitude. So it is not surprising 
that perturbing the first order solutions fails to capture the angular dependence of the second 
order solutions. 

Physically what the absence of the first order monopole from the second order Boltzmann 
equations means is that if we have a uniform radiation field then scattering by a stationary 
inhomogeneous distribution of electrons does not introduce additional inhomogeneities in 
the radiation field (in the elastic Thomson scattering limit). The dipole seen in the electron 
rest frame contributes to the additional inhomogeneities in the radiation field but it is 
small during recombination. Our analysis justifies neglecting the second order monopole, 
dipole and quadrupole, as we did in KW09. In particular, we conclude, as in KW09, the 
confusion with the primordial non-Gaussianity of local type resulting from inhomogeneous 
recombination is \Jnl\ ~ 1 and thus not important for the Planck satellite mission [3j which 



is predicted to achieve an accuracy of A/jv L ~5 [30j,|3l|. The S/N for the detection of this 
bispectrum by an ideal full sky experiment using temperature data alone is ~ 1. However 
perturbations in the electron number density will also have an effect on CMB polarization. 
If this effect is of a magnitude comparable or larger than the effect on temperature, a post- 
Planck, high-resolution, all-sky mission measuring the CMB temperature and polarization 
anisotropies may see the imprint of inhomogeneous recombination in the CMB bispectrum 
at few sigma level. 
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Appendix A: Contribution from Og^ and Vm 

We can write the formal solution for B^^k, n, t]q 
6( 2 )(k,n 

We will first include only the first term in Equation 1291 in the source S^(k, n, n). The 
calculation for other terms is similar. 

The angular averaged bispectrum is defined as sum over the m's of bispectrum times a 
Wigner 3jm symbol. 



,r] )= I d?7e^-^ k -V r S( 2 )(k,n,77). (Al) 
o 



B 



E 



1 <-2 t 3 



ITLi TTl2 7^3 



'1 *>2 «3 



E I" 1 " °\ ( a ^mi( x ^o)a^ m2 (x,T]o)afJ m3 (x,T] )) + 2 permutations 

m, 1 m 2 m 3 \ nil 1^2 ^3 

(A2) 
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,( 2 ) 
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where is the Fourier transform of 0^ and is calculated from first order multipole 
moments 09^. 



l(2) 



(1) 



Proceeding as in KW09 we get for the bispectrum from the first term in Equation [29] 



d 3 fc 



V mmi W \ J y Q (27T)3(27T)3(27r)3 l ^ 



P(k x )P(k 2 )(Anfl 2 f n d V '^k 3 c s ( V ') sin [fc 3 (r,fa) - r. (?/))] ji 3 [h( V - Vo )} 
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**„*&)l&(-k)<M*i,nOe^ + k 2 + fc s ) 



+ 5 permutations 



(A3) 



We have ignored IT* 1 ) in Sj$ to simplify equations, including it at the end of the calculation 
is trivial. We now use the Dirac delta distribution to integrate over k 3 . 



mim2m3 
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(An) 2 I " drjg(rj) 
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+ k 2 |c s (r/) sin [|k x + k 2 |(r s (77) - r s (r]'))\j h [|kj + k 2 1 (r? - t/o)] 

^ mi (^)r/ 2m2 (^ 

+ 5 permutations (A4) 



To proceed further we will need the following addition theorem for spherical waves 32] 



z L {\k! + k 2 |r)F LM (k 1 Tk 2 ) = J2 i ei+e2 ~ L (-l) M ^n(2L + l)(2£ 1 + 1)(2£ 2 + ^(M^M 

( a a 



li h L 




L 



^imi(^l)^ 2 m2 (fc 2 ), 



y mi m 2 —M 

where zg is any of the spherical Bessel function and the sum is over all allowed values of 
£i,£ 2 ,mi,m 2 . The above equation is valid for arbitrary values of k\ and k 2 if zt = je, the 
spherical Bessel function of first kind. If z# = t/£, the spherical Bessel function of second 
kind, then Equation IA5I is valid for k\ < k 2 (and for k 2 < k\ after interchanging k\ and k 2 ). 

We now use IA5I for the product ji s Y^ m . We also write sin [|ki + k 2 |(r s (?7) — r s (i]'))] = 
[|ki + k 2 |(r s (?7) — r s (r)')))jo [|ki + k 2 | (r s (r)) — r B {rf))\ and use Equation IA5I again. We also 



(A5) 
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use 



|ki 



k\ + k\ + ^Afc^Y^kOy^ka) 



(A6) 



The angular integrals over ki and \s.2 can now be done. Right hand side of Equation IA6I 
consists of two terms: k\ + k\ has no angular dependence while the rest of the right hand 
side depends on the angles ki and k2. For simplicity we will show the calculation for only 
k\ + k\ part. The calculation for the other part is similar but since we have extra factors 
of spherical harmonics we will get extra Wigner 3jm symbols on integration over angles 
summing over which will require few extra steps. 
The result for k\ + k\ part is 

(4tt) 3 
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(A7) 



Summing over the m's we get (32] 
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(A8) 



where the matrices in the last line are the 6j symbols. All the m dependence of the bispec- 
trum is in the above expression. Therefore to calculate the angular averaged bispectrum we 
need only consider the above expression for averaging over mi,m 2 ,m 3 . The result of doing 
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this averaging is 
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(A9) 



The calculation for the other term in IA6I is similar and it also results in the Kronecker delta 
symbol 6 20 = 0. 

The second term in Equation [29] involves cosine which can be written in terms of the 
spherical Bessel function of the second kind, 7/0 • We therefore need to break the integral 
over (k\, k 2 ) in two parts , k\ > k 2 and k\ < k 2 in order to apply the addition theorem. Both 
the terms will give a zero contribution to the angular averaged bispectrum (with Sio in the 
final result due to Ylo in this term), which is easily shown by a calculation similar to above. 
The boundary k\ = k 2 will also give zero contribution to the (ki,k 2 ) integral because the 
integrand is finite. 

(2) 

Thus we have shown that the contribution from 0q O to the angular averaged bispectrum 
vanishes. A similar calculation for the Vff shows that the contribution from the terms 
involving Ss2 in Equation |34] also gives zero contribution to the angular averaged bispectrum. 
In general Qlm ~ 8e®i Ytm gives non-zero contribution to the angular averaged bispectrum 
if and only if L — £ and M = m because of the orthogonality of spherical harmonics of 
different orders. 
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Appendix B: Integral Equation for second order monopole 



An alternative to solving the Boltzmann hierarchy for the second order monopole is to 



solve an integral equation 
equation is 



33 



341 ] . The line of sight solution for second order Boltzmann 
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where fi represents a general first order term multiplying S e . We can integrate over direction 
n to get an integral equation for the monopole 

eg ) (r/,k) = e^) fdrfg(rf) [(4tt) 3 / 2 !^-^(k!) £ j t , [k(rf - r,)] h„{k - k', rf) 
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We can now write down the contribution of O 00 to the bispectrum 
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Here we have used the integral equation for 0qq (Equation IB2j) to get an equation for 
S™j£l* m3 . The last term involving Q^m" wm §i ve a sman contribution (~ 10%) because of 
the factor of 1/10 and can be neglected. For vffl we can use the approximate tight coupling 
solution, the last term in Equation [3H in which case it can be absorbed into fan for I" — 1. 
We can similarly absorb the last term also if we choose not to neglect it. If we did not have a 
factor of jo multiplying the second order monopole term in last but third line, we would have 
an integral equation for S™],^™ 3 - We can however make progress by using the approximate 
solution for the second order monopole Equation [291 Then a calculation similar to Appendix 
lAl shows that the contribution of this term to the reduced bispectrum is exactly zero, so 
this term can be dropped. For the other terms we proceed as in KW09 and Appendix [A] 
We break the four point correlation function of first order terms into two point correlation 
functions using Wick's theorem. We can then perform all the angular integrals and two of 
the radial integrals using the properties of Dirac delta distribution, spherical harmonics and 
Wigner 3jm and 6j symbols. The result is 
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(B4) 

Note that this solution is approximate but does not assume tight coupling, despite the 
fact that we used the tight coupling solutions Equations [29l [33] and [34] as a trial solution. 
Equation IB4I is the result of iterating the integral equation once and will therefore contain 
corrections beyond the tight coupling approximation. In particular this solution takes into 
account all the terms in the full Boltzmann hierarchy, Equations [rTTfTQl The dominant 
contribution would come from around the last scattering surface, that is when rj' — 77 ~ 0. 
In that case the corresponding spherical Bessel functions would be close to zero unless the 
order of the spherical Bessel function is zero. Thus we would expect that most contribution 
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comes from terms with £'{ = = 0. The last Wigner 3jm symbol then forces I" = 0. But 
fen =0 = since the first order monopole cancels out making S™}^™ 3 vanish. This is the 
result that we found for the approximate solution of the second order Boltzmann equations 
also. For £' 2 ' ^ we also note that the arguments of the first two spherical Bessel functions 
differ from the arguments of the last two spherical Bessel functions by a factor of ~ 100. 
But for the squeezed triangles we would expect either t x or £ 2 to be small making £[ ~ £' 2 ' 
or £' 2 ~ £'2 due to triangle conditions in Wigner 3jm symbols. Thus we have a product of 
the spherical Bessel functions of similar orders but with arguments differing by a factor of 
hundred. This product will be negligibly small, since if one of the spherical Bessel function 
is near the peak the other would be negligibly small or oscillating very fast giving a small 
residual after integration. Thus the contribution from the second order monopole can be 
safely neglected for the case of inhomogeneous recombination. This argument also applies 
to all other terms in the second order Boltzmann equation which are a product of monopole 
type term and higher order multipoles. 
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